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Abstract 

We derive a general expression for the Casimir energy correspond- 
ing to two flat parallel mirrors in d + 1 dimensions, described by non- 
local interaction potentials. For a real scalar field, the interaction with 
the mirrors is implemented by a term which is a quadratic form in the 
field, with a nonlocal kernel. The resulting expression for the energy is 
a function of the parameters that define the nonlocal kernel. We show 
that the general expression has the correct limit in the zero width case, 
and also present the exact solution for a particular case. 

The increasing interest in the Casimir effect [I] is a natural consequence of the 
availability of new precision experiments, vv^hich pose an important pressure 
to continuously refine and improve the existing calculations. This evolution 
manifests itself at different levels; one of them amounts to coping with sit- 
uations where the geometry of the mirrors is more complex, albeit with an 
idealized description of their material properties, i.e., they are regarded as 
mathematical surfaces occupied by perfect conductors. Another level, which 
has recently received much attention, is the use of a more accurate descrip- 
tion of the mirrors, including corrections that represent their departure from 
exactly conducting surfaces: rugosity, finite temperature and conductivity, 
as well as finite width. The latter is usually dealt with by the introduction of 
a 'space dependent mass term' whereby the scalar field becomes very massive 
at the locii of the mirrors [2J. For the case of a real scalar field ip in d + 1 
spacetime dimensions, and a single flat mirror centered at Xd = 0, the local 
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Euclidean action, Siocai, for this kind of term is: 
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Siocaii^) = 2 y ^ y dXdV,{Xd)[^{x\\,Xd)] , (1) 

where x\\ denotes the time (xq) as well as the d — 1 spatial coordinates parallel 
to the mirror (which we shall denote by xy). The local potential Ve{xd) 
is a positive function concentrated around 0, with a width of size e. The 
perfect mirror case is approached when that size tends to zero and its strength 
becomes infinite; namely, ^(^^d) Q^i^d)- 9 ^ oo. This limit is a delicate 
step, since it usually introduces divergences [3j that may be difficult to deal 
with in a setup that used idealized boundary conditions only. 

Local interaction terms have been extended to include a non trivial depen- 
dence on the parallel coordinates [HE]. Translation invariance along them, 
necessarily implies a spatial nonlocality: 

S'locaii^) ^ ^ J ^'^^W J ^'^^'w J dxd^pixii,Xd)V,{xii- x\;Xd)^pix\,Xd) , (2) 

but a Fourier transformation in xj| yields a local expression in the mixed 
momentum (fc||) and coordinate (xd) representation: 
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S'locaAv) = 2 j (2^ y ^^df*{k,Xd)V,{k\i-Xd)f{k\\,Xd), (3) 

where the tildes denote the Fourier transformed of the corresponding object. 

Note that the resulting interaction term, S'i^^^i{ip), is still assumed to be 
local in Xd, and Ve{k\\;Xd) is concentrated around Xd = 0, on a region of 
size ~ e. However, except for the case of a zero-width mirror, a potential 
which is local in Xd can only be an approximate description of the interaction 
with a real material. Indeed, as explained in [6], one should in general use 
interactions that also include 'spatial dispersion' in Xd, i.e., nonlocality in 
the normal coordinates. As also shown in [6j, in spite of the nonlocality 
of the interaction, one may nevertheless use a Lifshitz formula [9] for the 
Casimir energy, since that formula depends on the reflection coefficients at 
the media surfaces, and they may be deflned also for nonlocal media, under 
quite general assumptions. Interesting results have also been obtained using 
the boundary state formalism [7j. 

We note in pass that the relevance of nonlocal media to Casimir-like 
effects has also been appreciated in other contexts, like String Theory [8]. 

In spite of the validity and usefulness of Lifshitz formula for nonlocal 
media, we believe that it would be important to have an alternative expres- 
sion for the Casimir energy, where its dependence on the details defining the 
nonlocal media were more explicit. 
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One might also be interested in situations where the material is not ex- 
actly confined to the region between two surfaces, but rather is concentrated 
on a region, with a non zero (albeit rapidly vanishing) density outside of that 
region. This is a situation where the use of reflection coefficients, albeit still 
possible, becomes nevertheless problematic. 

To confront those problems, we first define the setup: for a mirror centered 
at an arbitrary position Xd = b, we shall use an interaction term Sf'^ given 
by: 

X i^*{k\\,Xd)Ve{k\\]Xd-h,x' d-h)^{k\\,x d) ■ (4) 

The kernel V^{k\\]Xd,x' d) is not invariant under translations in the normal 
direction {xd —>■ Xd + h, x'^ ^ x'^ + h); besides, we still assume it to be 
concentrated around Xd = and = 0. This concentration may be used 
to write a convenient expansion for the kernel, based on the introduction 
of [^n\xd)]^, an orthonormal basis of functions of the normal coordinate, 
obeying the boundary conditions that follow from the microscopic model 0. 

Then, without any loss of generality, the nonlocal kernel will be expanded 
as follows: 

V,{kfXd,x'd) = J2 Cmn(fc||,e)^^Ha^d) ^i'^*ix'd) , (5) 

m,n 

where Cmn{k\\,e) = C*j^{k\\,e), from the reality of the action. 

Being a nonlocal term, one has to rephrase the properties it should have 
to behave as a (generalized) mass term, in the sense that it favours the 
vanishing of the field around the region where it is different from zero. It 
is clear then that the quadratic form ([5|) has to be definite positive (in the 
space generated by the basis); this amounts to a non-trivial relation for the 
Cmn matrix. 

Let us illustrate the previous construction with two examples. Firstly, 
we consider a model where Sf^ emerges from the linear coupling of (/9 to a 
microscopic real scalar field ^(x), which is confined to the region jx^ — 6| < e/2 
(xy: arbitrary), satisfying Dirichlet boundary conditions at x^ = 6±e/2. It is 
sufficient to deal with 6 = 0, since the general case is obtained by a translation 
of the kernel. Following a generalization of the approach of [lO], we see that, 

^For example, they could satisfy Dirichlet boundary conditions at Xd = ±e/2. But they 
are not necessarily of compact support; they could for example be exponentially decaying 
functions with a typical dispersion ~ e 
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in the functional formalism, Sf^ may be written as follows: 

r Tli: ^-Smi^) + ig f d'^x\\ J+'/j ^.^d .x^) 

where (7 is a coupling constant and 5*^ is the action for the microscopic field. 
The matter field ^ may have a self-interaction, controlled by an independent 
coupling constant (implicit in Sm)- 

To proceed, we denote by W{J) the generating functional of connected 
correlation functions of ^, related to 2{J), the one for the full correlation 
functions: 

hy W = InZ. The current J is confined to the same region as ^, but it has 
free boundary conditions. We then have that Sj{(p) = —W[igip(x)]. On the 
other hand, since only the quadratic part in (p will be retained □, 

Sf{^) = -W[zgv{x)] 

- ^9^ J d'^+^x J d''+^x'ip{x)W^^\x,x')ip{x') , (8) 

where W^"^^ is the connected 2-point function. 

Now, assuming that 5*^ is translation invariant along xy, we immediately 
identify the nonlocal kernel: 

= g^W^^\k^^-x,,x',) . (9) 
Besides, since C, satisfies Dirichlet boundary conditions, we have: 

W^'\kf±e/2,x',) = W^'\kfx,,±e/2) = 0. (10) 
Then we have the expansion: 

V}'\kfX,,x',) = J2^l^\^d)CUh^^)i^ln^*i^<i) (11) 
m,n 

where the orthonormal functions are given by: 

-(e). . _ ^/2 r sin(^) if n = 2k, {k = l,2,...) 



^The media are assumed to be linear. 
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The precise form of Cmn{k\\, e) depends on the action Sm- If it is a free 
action we have the diagonal expression: 

Cmn(fc||,e) = ^nny _^ ^2 _^ ^2 ' ^^^"^ 

where /i is the mass of the microscopic field. 

As an alternative example, we consider the case of a charged field ^, ^, 
coupled quadratically to the real field (p. In this case, the analogue expression 
to ^ would be: 

e-^f (V) = J ^^^^^ 

JV^V^e-s^i^'O ' ^^^^ 

where, to simplify the treatment, we assume Sm to be quadratic: 

Sm{^,0 = [ d% r ' dx^{didi + /i^a] • (15) 

J J-~e/2 

The integral can be formally performed, but the result is a functional deter- 
minant. Since we use a quadratic approximation for Sf \ we only need it up 
to that order in (p. The corresponding contribution is just a 1-loop diagram 
with two legs. Translation invariance along the parallel coordinates suggest 
the use of a mixed Fourier representation: 

G{p\\ ; Xd, x'JG(p|| + k\\ ; x'^, Xd) (16) 

where G{j>\\] Xd, x'^ is the microscopic field propagator in the mixed represen- 
tation. Rather than evaluating the actual form of the nonlocal term for this 
model, we just want to show that the boundary conditions for the nonlocal 
kernel are determined from the ones we impose on the microscopic field: as- 
suming, for example, that this field satisfies Dirichlet boundary conditions 
at the boundaries of the mirror, from (fT6l) we derive for the same kind of 
condition. Namely, the kernel vanishes when Xd = ±e/2 or x'^ = ±e/2. Thus, 
also in this case the model produces a nonlocal kernel with the structure of 
(fTTl) (the same basis), with different coefficients C^n- 

Before evaluating the energy, let us briefly examine the boundary condi- 
tions that follow from the nonlocal term, in a concrete case. To that end, we 
consider the real-time version of the equations of motion for a free massless 
scalar held coupled to a nonlocal potential centered at Xd = 0. Assuming, 
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for the sake of simplicity, d = 1, and Cmn = Cmn{^) (independent of cj), the 
equation of motion becomes: 

(17) 

Fourier transforming in time, 

{-dl - u;^)<^(u;, x'J = - / dx'^ V,{xi, x\) (p{uj, x\) . (18) 



Then we multiply both sides of the equation above by il)m*{xi) and integrate 
over xi, to obtain: 

{i^^^\{-dl-uj')\m) = -$^C_(6)(^(^)|(^(c.)), (19) 

n 

where Dirac's bracket notation denotes the scalar product of functions of Xd- 
Let us assume, for the sake of simplicity, that the proper basis is (fT2l) . Then 
{il)m\{—d1)\!f>) = {^)'^{ipm\^) and, as a consequence: 

J2 Cmnie) an = [u' - (^)'] a^, (20) 

n 

where a„ = {ijjn^ \(p{k\\)) . Thus, 



J2 Cmn{e)alan = J] - i^Y] . (21) 



This means that, to have a solution, the a„ coefficients diagonalize C, and, 
since the quadratic form on the left hand side is positive, we am vanishes 
whenever tu'^ < {^)'^. 

This means, in particular, that for < (7)^, all the coefficients vanish: 
the field vanishes when \xi\ < |. Of course, things are different if, for exam- 
ple: (-)^ < u"^ < (— )^, then only the first coefficient may be different from 
0. 

An interesting case is that of a Cmn which is a finite matrix: one that 
vanishes for m > of , say. An extreme case is A^ = 1: there are then only 
two regimes, depending on whether u"^ is bigger or smaller than (^)^. In the 

former case, (p in orthogonal to This implies that has at least one node 
in the [— f 7 f] interval. This is the manifestation of a Dirichlet-like boundary 
condition in this context, which of course will only hold true up to certain 
values of u"^. For bigger values, the previous condition is relaxed and the 
mirror is transparent. 
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Let us now evaluate the Casimir energy, discarding terms that are inde- 
pendent of the distance between mirrors (and do not contribute to the force). 
For two mirrors, one at = and the other at Xd = a, the total action S is 
S{ip) = So{ip) + Si{^) where So = \ j d'^+^x d^ipd^^, and Si{^) = + 

Since S* is a quadratic in the fields, it is immediate to find an expression 
for the vacuum energy Eq, in terms of the determinant of the corresponding 
kernel defining the quadratic form. Since we have translation invariance along 
X||, we use the energy per unit area, Sq, and take advantage of the Fourier 
transformation to obtain: 

1 f d%\ 



where K, is an operator acting on functions of Xd, whose matrix elements are: 

K,{xd,x'^) = }Co{xd,x'^) + V}°\kfXd,x'd) 

+ V}''\kf,Xd^x'd) . (23) 

where }Co{xd,x'^) = {—d1 + /cy)5(xd — x'^. The trace operation, denoted by 
'Tr' refers to the trace in the space of functions depending on Xd- 

The expression above contains three contributions which, to calculate the 
Casimir force between the two mirrors, are irrelevant. One of them, 8^°''^, 
corresponds to the vacuum energy density in the absence of mirrors 

The other two, denoted by £^0°^ and £'^\ are the mirrors' self-energies (and 
therefore we shall discard them). They have the form: 

= 5 + "'■"')• <^^> 

where 6 = 0, a, and W^'^^ denotes an operator, acting on the same space as 
above, and defined by: ^ 

= /C-i k(^) . (26) 

^(0)^ y{a) l^j^^g 

non trivial matrix elements in smaller spaces, namely, the 
ones generated by the basis functions sitting at each mirror. This fact can 
be used to show that the trace operation above may be taken, for each term, 
using only the respective basis at = and Xd = a (the trace operation 
over the complement vanishes). 



7 



Using matrix elements defined with the functions ipn^ and 0n "* , such that 
(pl^^ {xd) = i>n^ {xd - a) : 

= {<l>^:.^\IC^'V^'^^\<f>^:^) . (27) 

I denotes the identity operator, so that: Imn = ^mn- A simple shift of 
variables leads to Sq^^ = Sq°'\ as it should be. 

After extracting the previous three contributions, we obtain a subtracted 
energy density, Sq, which by some straightforward algebra may be put in the 
form: 

where O is an operator whose matrix elements may be given in terms of the 

C^mn = ^ ^ ^rnp^pq^qr ^rn y (29) 



i/jn '' basis, as follows: 



p,q,r 



where 



U^l = (0^)|[/Co + K^'^)]"^!^^) = U^l ^ Umn. (30) 
Taking into account the previous relations, 

which is a nonlocal version, for flat, identical mirrors, of the equation derived 
in [m (see also [T2]l 

Note that Umn can be written more explicitly, by performing an expansion 
in powers of Deflning: A^^ = {^'ti\\Ki<:^~^\\\}'n) and r^n = (■?/'-m^| [/Co]~"^|0n^), 
we see that: 

7-r(0) _ Y — f\ r V 

^ mn ''nn '-^mp^pq^ qn 

-\- ^ff^pCpq^q^-C^gV -|- . . . (32) 

(we used Einstein's summation convention). Then we see that: 

U = (/ + AC)"^r . (33) 
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Let us check that the expressions above do yield the proper answer when the 
hmit corresponding to the case of perfect mirrors is taken. This is, in the 
present formalism, tantamount to: 

V;(°)(fc||;x,,x',) = \6{x,)6{x',) , (34) 

(sharp boundary conditions) and then A — * oo (strong boundary conditions). 
One immediately gets, from (i34l ). the matrix elements: 

C^n = A^^)(0)^(^>(0). (35) 

Inserting this into (i29l ). we obtain: 



1 r d%\ 



£-0 = — I -, — r-j In 



which in the strong limit, and for d = 3 yields: 

- d (4^"('-^ " 

7r2 



1440 a3 



(37) 



which is the proper result. 

Finally we consider a truly nonlocal example: 



Cmn{k\\,e) = 6,noSnoKh) ^ (38) 

where A(A;||) > 0. It corresponds to a 1 x 1 matrix C in the space generated 
by the basis functions. In this case, we find that: 

Omn = [Kk\\)Y UooUmoSnO ■ (39) 

For this particular case, the O matrix has rank 1, so only one of its eigenvalues 
is different from 0, and the trace of the log may be evaluated exactly. The 
Casimir energy becomes: 



1 r d^k 



^° = 27 (2^ 



>.Hk,){U„y- . (40) 



We may produce a more explicit expression for the matrix elements of U : 

U^^= [l + A(A;j,)Aoo]"'roo, (41) 
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where: 

dkd {{-ipo^lkd) 



Aoo = / ^ (42) 



2n k\ 



\kd) is the plane wave ket). 
The remaining object, Too is: 



00 



2vr kl + kl ■ 



Equipped with the previous expressions, we may calculate the Casimir energy 
for a particular basis element. An interesting example is the choice of the 
exponentially localized function 



= ' (44) 

since it allows one to evaluate the integrals above exactly. The result for £q, 
in d spatial dimensions, may be written in terms of an integral involving a 
function G 



(45) 



r (f ) Wl'^a'^ Jo 

depending on dimensionless parameters. It is given explicily by: 

Gip; X, /) = ^-''-M^P-p+^'-Mxp+^J 

^ + 2erfc(xj>) 

where erfc is the complementary error function. 

It is immediate to check that the Casimir energy for the perfect mirror 
case is reproduced when ^ — > and A -^^ oo. On the other hand, when that 
limit is taken for a finite A, the result has the same form as in the local case, 
but with a (finite) renormalization for A. Indeed, if Xiocai denotes the coupling 
constant in the local (5-potential case, the energies agree for A = 

For finite values of -, we have the interesting phenomenon that the cor- 
rections are not analytical in that variable. Indeed, one can see that the 
corrections to the zero-width case are proportional to a factor e~^^^ . 

To make the comparison with the local case more explicit we present, in 
Figure 1, the plots of the Casimir energies corresponding to the local and 
nonlocal cases, for d = 1. The local potential is chosen so that it agrees with 
the nonlocal one when e ^ {Xiocai = 8^/^TX)■. 
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Figure 1: Casimir energies corresponding to the nonlocal (continuous line) 
and local (dashed line) cases, as a function of 6 = | for eA = 



A remarkable fact, that can be observed in the plot, is that the Casimir 
energy for the nonlocal case becomes finite when the distance between the 
mirrors tends to zero. This is a manifestation of the fact that the nonlocality 
softens the UV behaviour of the system. Yet another consequence of the 
same effect is that the integral over k\\ for the energy of a single mirror, ^q'''' 
(see Eq. (l25ll). has a better UV behaviour that its local counterpart. In 
particular, for ci = 1, a simple calculation shows that it becomes: 



£^^'^ = — I dk In 



1 



'0 



27r 







1+^ e'^''\Tfci2ek) 
k 



(47) 



which is convergent for any e > (we recall that its local counterpart is 
logarithmically divergent [3]). 

We conclude by noting that, as shown in the examples above, nonlocal 
potentials can be used to impose boundary conditions in finite size mirrors, 
and they becomes automatically frequency dependent. Also, in spite of their 
seemingly complex structure, a general expression for the energy may be 
derived, which contains some new interesting features: the non-analytic be- 
haviour of its small-width expansion and a softer UV behaviour. 

In spite of the above, the perfect mirror limit is still properly reproduced. 
Besides, when the distance between mirrors is of the order of e, the Casimir 
force vanishes, rather than becoming infinite, as it happens in the local case. 
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